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Abstract. We present a general model for an atomic memory using ultra-short pulses of light, which allows 
both spatial and temporal multimode storage. The process involves the storage of a faint quantum light 
pulse into the spin coherence of the ground state of yl-type 3-level atoms, in the presence of a strong driving 
pulse. Our model gives a full description of the evolution of the field and of the atomic coherence in space 
and time throughout the writing and the read-out processes. It is valid for any frequency detuning, from 
the resonant case to the Raman case, and allows a detailed optimization of the memory efficiency. 
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1 Introduction 

For quantum telecommunications and for quantum infor- 
mation processing, memory registers able to store quan- 
tum information without measuring it are essential de- 
vices. A quantum memory relies on an efficient coupling 
between light and matter, allowing reversible mapping of 
quantum photonic information in and out of the material 
system. In the past years, several protocols have been de- 
veloped theoretically and experimentally [HIS]. Storage and 
retrieval of some of the basic states of light for quantum 
communication such as a polarization q-bit [3 , squeezed 
light [HIS] and entangled photons [6;,7| or faint coherent 
pulses at the level of one to few photons have been re- 
alized [SlI^ fTUlimiT^ ■ Recent experiments have achieved 
large efficiencies [TUllllj . 

However, the processing speed and the available band- 
width of the memories remain a challenge for quantum 
memories. The first quantum memory registers proposed 
more than a decade ago [13| involve the transfer of quan- 
tum information from light to atoms (writing) and back 
from atoms to light (retrieval), using electromagnetically 
induced transparency (EIT) in atomic three-level transi- 
tions, and this process implies a limited bandwidth. The 
storage protocol relies on a strong control field, generat- 
ing EIT for the weak field that carries the quantum sig- 
nal to be stored. The group velocity for the signal field 
is strongly reduced and the signal pulse is compressed by 
several orders of magnitude. A signal pulse can thus be 
contained inside the atomic medium, and before it propa- 
gates outside the medium, the control is switched off. The 
quantum variables of the signal field are then converted 



from a purely photonic state to a collective spin coherence. 
For read-out, the control field is turned on again and the 
medium emits a weak pulse, carrying the quantum infor- 
mation contained in the original pulse. While in principle 
this allows direct mapping of the quantum state of light 
into long lived coherences in the atomic ground state, the 
bandwidth of the stored signal is strongly limited by the 
transparency window associated to EIT. 

Various methods have been proposed to achieve broad- 
band memories and escape the limitations linked to EIT. 
Very interesting protocols are based on the implementa- 
tion of controlled broadening, such as CRIB (controlled 
reversible inhomogeneous broadening) |14| . using photon 
echo-type reversal ^15,, and AFC (atomic frequency comb) 
where an absorbing comb structure is created in the medium 
[IB] . These techniques have been successfully applied for 
echo-type light storage in rare-earth doped crystals [TT] 
and atomic vapours [18) . These methods allow broad band- 
width but they imply writing times which are still rather 
long (of the order of microseconds) . In the spatial domain, 
an interesting phenomena, quantum holography, has been 
proposed to implement BD-memories |T^I20j . 

An alternative method is based on the use of a broad- 
band, ultrafast control and signal field pulses for the writ- 
ing process. It relies on two different atomic transitions 
sharing the same excited state. The control field and the 
signal field contribute to a two-photon process coupling 
two ground states. In this case the signal pulse is con- 
verted into an atomic coherence between ground and ex- 
cited states and then into a ground state coherence by 
the control pulse. However, since the interaction times are 
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Fig. 1. Three level atomic system interacting with driving field 
n and signal field a. 



very short, it does not allow for the buildup of EIT. This 
method has been proposed pill^ and demonstrated ex- 
perimentally in far ofF-resonance conditions [251IT^ . 

In this paper, we present a detailed theoretical model 
for an ultrafast memory without adiabatic approximation 
and valid for arbitrary frequency detuning. We show that 
this protocol holds the promise for a quantum memory 
with high efficiency, fast operation and broad bandwidth 
together with spatial multimode capacity. The problem 
of the achievable efficiency in this case was treated in 
Refs. |2ip22j using a numerical optimization procedure 
based on the search of the optimal pulse shape for the 
signal or driving field in the limit of adiabatic elimination 
of the excited state. In Ref. [23] the shaping of the driving 
pulse is based on the analysis of the Lagrange function, 
avoiding the adiabatic approximation in the optimization 
procedure. A very good efficiency can be obtained even 
for short pulse durations. In Ref. |25| a different optimiza- 
tion technique based on the minimization of the losses has 
been used in the resonant case without adiabatic approxi- 
mation, in the limit of very short pulses (shorter than the 
excited state decay time). In the present work we extend 
this technique to the case of arbitrary frequency detuning. 
We show that the memory efficiency can be very good even 
for large detunings as long as the experimental parameters 
are properly optimized. Moreover, we explore the transi- 
tion region between resonant and adiabatic regimes and 
we demonstrate that our technique allows finding optimal 
parameters for storage. 

The article is organized as follows. In Section 2 we 
present the model system and we write the main equations 
ruling it. In Section 3 we give the method for solving the 
equations for the writing and read-out processes in the 
semi-classical limit. In Section 4, we study the evolution 
of the atomic coherence and of the signal field during the 
writing process and we calculate the losses. In Section 5, 
we study the read-out process and the efficiency of the 
memory as a whole. 



2 Model system 

In this paper, we consider an ensemble of three-level atoms 
in a yl-configuration (Fig.[T]) that will be used to store tem- 
poral and spatial multimode quantum fields. The atoms 



interact with two electromagnetic fields, a signal field 
and a driving field i?^, that connect the two atomic ground 
states to the excited state. The driving field is a strong, 
classical field propagating as plane wave, while the signal 
field is a weak quantum field with a transverse structure. 

The signal and driving field are very short pulses that 
are assumed to be much shorter than the excited state life- 
time 7~^, so that we can neglect the spontaneous emission 
during the writing process. 

In the dipole approximation the light matter interac- 
tion Hamiltonian is given by 



V = - 



3 



dj{t)E{t,r,), 



(1) 



Here dj (t) is the electric dipole operator of the j-th atom 
located at rj . In the paraxial and quasi-resonant approx- 
imations the Hamiltonian can be rewritten in the form 

JJJdzd^p ihg d{z,p,t)a3i{z,p,t) e'^'^"^^'^^^ 
'f2{t)<^32{z,p,t) e'k,z-zAt 



Ah 



-f2*{t)a23iz,p,t) e 



-ikdZ + lAt 



(2) 



Here kg and kd are the wave vectors of the signal and 
driving fields, and p = p{x,y). The one-photon detunings 
of the signal and driving fields are assumed to be equal, 
and equal to A so that the two-photon resonance between 
levels 1 and 2 is fulfilled 



(3) 



The spatial coordinates z and p describe the longitudinal 
and transverse propagation of the signal field. The nor- 
malized amplitude of the signal field a(z,p,t) is written 

as 



E,{r,t) = -i 



2eoc 



a{z, p,t) + h.c, (4) 



where a{z,p,t) is the annihilation operator for the signal 
field. Under propagation in free space we have the follow- 
ing commutation relations [26' 

[d{z, p, t),d\z, p', t')] =S'{p- p')S{t - t'), (5) 
[a{z,p,t),d\z',p',t)] = 

The amplitude d{z,p,t) is normalized so that the mean 
value {a^ {z, p,t)a{z, p,t)) is the photon number per sec- 
ond per unit area. The symbol A± represents the trans- 
verse Laplacian with respect to p 



A, 



(7) 
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We define the intensity of the driving field through the 
Rabi frequency Q. For the sake of simplicity we consider 
this value as real fl — Q* . The driving field is a classical 
plane monochromatic wave propagating along the z-axis. 
The coupling constant between atom and signal field is 



According to equations ( 10 ) the collective atomic op- 



/ \ 1/2 



(8) 



where d^i is the electric dipole element on the transition 
|1) — 1 3). We define the collective coherences and popu- 
lation as sums of over all the atoms 

a,u{r,t)=Y.^l,{t)5^{r-r,), (9) 
iV,(r,i) = ^a^,W<53(r-r,-)- (10) 



These quantities fulfill the commutation relations 

[(j,k{r,t),a,k{r',t)] = \N,{r,t)~Nk{r,t)\ S^r~r'), 

(11) 

In this basis one can derive a full system of Heisenberg 
equations for the collective operators namely the field am- 
plitude a{r,t), the collective atomic coherences aik{r,t) 
and the collective atomic populations Ni{r, t). Taking into 
account the Hamiltonian ^ and the commutation rela- 
tions (|6|, ( [TT| ), the complete system of the equations reads 



d 



d 



dt '^dz 



ic 

Ik a 



A, 



dt 

dt 

dt 

dt 
d_ 
dt 
d^ 
dt 



0-13 = -i/4(7i3 + Qai2 

(7i2 = -i7(5-i3 - 5(2(732, 



(T32 = iAa32 - ^^(^3 - ^2) + 5a^o'i2, 
^1 = -90.^31 - ga^ ^13, 

N2 = - f2 ((732 - (723) , 

dt 



d_ 
dt 



(12) 

(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



To derive the equations (12)-(18l we have performed the 
substitutions 



ik^z — iAt ~ 



0-13 -> e"-"- ■ "(713, 

» , JkdZ — iAt±. 
0^23 e (723, 

» ^ ^i{k,i — k^)z " 
ai2 — J> e " (7i2. 



(19) 



We have omitted the terms related to spontaneous re- 
laxation 1 3) — >■ |1) since we have assumed that the relax- 
ation rate 7 is small enough for the spontaneous emis- 
sion to be negligible during the short time duration of the 
pulses. The equations can be written in a simplified way 
with a few approximations given below. 



erators have sharp spatial distributions due to the delta- 
localization of the atoms. However, due to the collective 
effect of atoms localized along a field trajectory, we can 
average Eqs. (12)-(18l over the positions of the atoms. 



We will also replace the operator Ni — N3 in Eq. ( |13[ ) 
by the number giving the mean atomic density TV. We will 
assume that in the beginning of the process most atoms 
are in state |1). During the memory processes (writing 
and read-out) the population of the state |1) stays close to 
its initial value, because the photon number in the signal 
pulse is much smaller than the initial atomic number 



In Eq. (14), we can neglect the second term on the 
right hand side since ga is much smaller than f2 (we have 
assumed \Q\'^ ^ g'^{a^a)). Furthermore, (732 ^ (7i3 be- 
cause the populations N2 and iV3 are much smaller than 
N. 

Now we can write a simplified system of partial dif- 
ferential equations describing the evolution of the system 



ld_ d_ 
c dt dz 



2,ka 



V_L a,{z, p, t) = ~g 0-13(0, P, 0>(20) 



_9 

m 

d_ 
dt 



^uiz, p, t) = -iAai3{z, p, t) 

+gNa{z,p,t) + n&i2{z,p,t), (21) 
&i2iz,p,t)^-n&i3{z,p,t). (22) 



In a previous publication |25) . similar equations were de- 
rived for the case of a resonant excitation {A — 0). Here 
and below we omit the averaging over the atomic localiza- 
tions. 

Let us renormalize the coherences (713 and ai2 

(23) 
(24) 

so that they obey to the bosonic commutation relations: 



(5-12(2, P, t)/y/N ^ b{z,p,t), 
^i3{z,p,t)/y/N = c{z,p,t) 



bir,t),b\r',t) =[cir,t),c^ir',t)]=S^r-r'). (25) 



Here again we have taken into account the fact that Ni — 
N2.3 N. 

In the Fourier domain relative to the transverse coor- 
dinates p the equations read 



a{z,t:q) = —gVN c{z,t;q), 



d_ 

dz 
d 

—c{z,t;q) = -iAc{z,t;q) 



d_ 
Ft 



(26) 



+gVN a{z,t;q) + [2b{z,t;q), (27) 



b{z,t;q) = ~f2c{z,t;q), 



(28) 
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where we have introduced the transverse wavevector q and 
we have made the changes 



6(z,t;q)^S(z,t;g)e-*9V(2M, 
c(z,t;g)^c(z,t;q)e-*«V(2fc.). 



(29) 
(30) 
(31) 



From the system (26)-(28) one can obtain a conservation 
equation 



da^a dWb 
+ 



dz 



dt 



1h 



0. 



(32) 



This equation means that the input photons of the weak 
quantum field are converted into excitations of the atomic 
coherences (Ti3 and (J12. The aim is to store the informa- 
tion carried by the signal field in the ground state coher- 
ence <7i2, so that an excitation of the state |3) is undesir- 
able. It is possible to reduce this loss channel by an appro- 
priate choice of the driving field. If the driving field power 
is high enough for the Rabi oscillation on the transition 
1 2) — >■ 1 3) to be more effective than a spontaneous emis- 
sion {{2 7) and if the pulse duration is short enough, 
so that the atoms undergoing a Rabi oscillation have no 
time to go back to the state |3), then the third term in 
Eq. (32) should be negligible. This will be studied in the 



optimization of the memory process. 

We neglect the time delay linked to the pulse propa- 
gation in the atomic medium. This means that, if we have 
long enough pulses, such that L/c T {L is the thick- 
ness of the medium and T is the pulse duration), we can 
neglect the time interval between the time at which the 
front part of the pulse enters the medium and the time at 
which the front part leaves it. Formally this means we can 
neglect the time derivative in Eq. ( 20 1 . For simplicity we 
will assume that the driving pulse has a rectangular time 
distribution (in the equations fl(t) — const for < i < T). 



3 Writing and read-out processes in the 
semi-classical limit 

The main aim of this paper is the determination of the 
memory efficiency. As it is well known the semiclassical 
description is sufficient for this and we can use the ini- 
tial conditions : &(0, z; q) — c(0, Z]q) = for writing and 
am(^; 9) = c(0, z]q) =0 for read-out. Here and below we 
omit the operator notation for the variables. The detailed 
resolution of the system of partial differential equations 

^13 



(26)- (28 1 can be found in Apps. 



Using the general solutions (|5ip-(53l one can obtain 
the semi-classical ones for the writing process for both 
the field amplitude a}^{t,z\q) and the atomic coherence 
h{t, Z] q) in the form 



a^{t,~z;q)^ / dPa,n{i',q)Gaa{i-i',S), (33) 
Jo 

b^il ~z- q)^~p dPa,n{i', q)Gab{i- i' , z), (34) 
Jo 



where we have introduced the dimensionless time t and 
longitudinal spatial coordinate z according to 



t = nt, 



Ty/ = f2T] 

n 



(35) 



where Ty/ ^ is the common duration of signal and driving 
pulses for writing and L is the thickness of the memory 
cell. We have defined an effective interaction coefficient p 
given by 



(36) 



The kernels Gaait, z) and Gab{t, z) are time convolutions 

Gaa{i,i)^ I di'f{i',i-r)r{t-i',z;~r), (37) 
Jo 

Gabii, z)= I dp Mi', ~z- r)f; [i - i', S; -r), (38) 
Jo 

where the functions / and /o are expressed via the n-th 
Bessel function of the first kind denoted by J„: 



/(t,z;r)= 5(t)-e-»^V*^l+^^~ 



XM^-^^^ .h[^Jil + r)2t]0it), (39) 



/o(t, S; r) = e~'^^e-'^^^ ^ J„ (^^ {1 + ryzi^ 0{i). 

(40) 

Here the frequency detuning is given by the dimension- 
less parameter r = A/{2S1) and 0{t) is the step function 
0{i) = 1 for < f < TV and equals zero otherwise. 

In order to estimate the efficiency of the memory, we 
calculate the field amplitude at the output of the cell. This 
amplitude is different in the case of forward and backward 
retrieval. The corresponding expressions read 



a%^{i,L;q) 



(41) 



- dt'a„i{i' ,q) j dzGabit,z)Gba{i',L-z), 
^ Jo Jo 



Gba{i,z) = Gab{i,z), 
and 



^back 



(t, L; q) 



(42) 



1 



Tw 



dt'a^n{i',q) / dzGabii,z)Gba{i' 



The last two formulas are correct only in the approxima- 
tion where diffraction is neglected. Indeed, we have not 
taken into consideration the diffraction phenomenon de- 
scribed by equation (29l-(31). As was demonstrated in 



[2 5) this effect does not introduce any significant correc- 
tions in the case of the forward retrieval but restricts the 
mode number for the backward retrieval. 



T. Golubeva et.al.: High speed yl-type memory with arbitrary frequency detuning 



5 



In the following we will use an optimization procedure 
based on the choice of the optimal relation between the 
thickness of the memory cell L and the duration Tw of the 
signal and driving pulses. We will compare this optimiza- 
tion approach with the approach developed in Ref. [22] 
based on the optimization of the driving pulse shape in 
the case of short pulses. 



4 Discussion of the writing process 

In this section, we study the writing process, i.e. the con- 
version of the signal field into atomic coherence. Let us 
start with a simple calculation of this process at the input 
of the memory cell. For this we solve Eq. ( 34 ) at i = 0. In 
this case the atomic polarization 0; qjcorresponding 
to the "written" information reads 



6*^(^,0) = 6^ (i,0;q)/(-2pa„(q)) 



- e 
ir 



SlU' 



-i{Wr^ + 1) 



iWr^ + 1) 



(43) 



For the sake of the simplicity we have taken ain{t,q) = 
const{t). The normalized atomic polarization ranges from 
to 1. In the limit of small or large detuning, we have a 
simple behaviour 



r<l: |6'^(t,0)r =sin* 



r>l: |6^(t,0)r =sin^ 



4r 



(44) 
(45) 



Fig. [2^ shows the Rabi oscillation for a normalized detun- 
ing r = while Fig. [2]: shows a periodical modulation with 
a period of inr for r — 10. When r increases, one can see in 
Figj2}3,c the beat between two close frequencies. For larger 
r (Figj2ji,e) the curves show a slow oscillation (that cor- 
responds to the term with frequency ^/r'^ + 1 — r in (43)), 



modulated by a fast oscillation (at frequency \/r'^ + 1 + r) , 
with a modulation depth decreasing with increasing r. 

Following the variation of |6'^(f, 0)p as a function of 
r allows to get a first view of the behaviour of the system 
when the detuning is varied. It can be seen from Fig.[2]that 
even for a small detuning, r = 0.1, a significant distortion 
of the coherence profile at the input of the medium takes 
place as compared to the resonant case. On the other hand 
for r = 2 the high frequency modulation of the slow oscil- 
lations is rather small and the excitation can be considered 
as close to the off-resonant case, where the well-known so- 
lutions in the Raman limit can be used. A detailed study of 
the coherence distribution for all z, as given below, will al- 
low to better characterize the interaction regime, between 
resonant and Raman. 

When shifting into the medium, the behavior of the co- 
herence is much more complicated since it looses its simple 
harmonic character. Figure [3] shows a displacement of the 



maximum of the coherence along the medium. Let us fol- 
low the dependence of l&'^i, z)P on z for a given value of 
i=f^=TT given in Fig. [t] 

Comparing the curves for different values of the de- 
tuning, we see that the value of coherence at z = sig- 
nificantly decreases with increasing detuning. From these 
curves it could be concluded that the larger r the less 
information is written in the cell. However, we will show 
that it is not true. Let us introduce the quantity n'^f^ that 
characterizes the proportion of signal photons converted 
into coherence during the writing process: 



1 1 



(46) 



Here, before the integral comes from the input pulse 

energy, and the factor 1/2 comes from the previously in- 
troduced dimensionless variables. The integral gives the 
normalized population in the medium with length L 
during the interaction time . Since the transition of an 
atom to the level |2 > in our model corresponds to the 
coherent scattering of a photon from the signal wave, this 
is also the number of signal photons that was recorded in 
the atomic coherence. We are interested in the percentage 
of input photons recorded in such a way. The calculated 
values of n'^-^-^ are presented Fig. [i^-c. We see that the 
number of recorded photons in the first two panels are al- 
most identical (there is even a small increase of n'^-f-^ for 
r = 0.5 as compared with the resonant value), while for 
r = 2 the value n'^^^ decreases by a factor of less than 2. 

We can also follow the dependence of n'^-l'^ on r for 
a given medium length and for different durations of the 
signal pulse. Fig. [5] shows that the proportion of recorded 
photons decreases when the detuning increases, but this 
decrease depends on the duration of the writing process. 
In particular, we see that the writing efficiency depends 
weakly on the detuning for = tt/4 but it stays quite 
low. With increasing pulse duration {T^ = tt) a plateau 
appears for detuning range to 0.7. This range is reduced 
for — 2tt, but a second plateau appears for r e [1, 1.8]. 
While Figs. [3j [4] and [5] give a detailed behaviour of the ef- 
ficiency of the writing process depending on the detuning, 
on the time duration of the pulse and on the length of the 
medium, it appears that the optimization of the efficiency 
is non trivial and requires a specific procedure. This will 
be studied in the next section. 



4.1 Estimation of the writing losses 

In Ref.[25j an algorithm of memory optimization, based 
on the minimization of leakage is described. Leakage is 
defined as 

£(f ^, L) = /o^" \a'^it,L)\^dt ^ ^^^^^^ ^^^^ 

It characterizes the proportion of signal photons going out 
of the cell during the writing time. Such an estimation of 
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the losses is justified when the leakage is the main origin 
of losses. However, as will be shown below, there is a range 
of in which the population of the upper energy level is 
large enough to cause significant losses. We will show that 
for a high-speed memory the three-level atomic system 
can not be reduced to a two-level scheme. Note that this 
situation is specific of the case of simultaneous interaction 
of the signal and control fields with matter, and does not 
happen in memory protocols based on an echo [15]. We 
have already introduced the value n'^-^-^ , which is the pro- 
portion of signal photons that have been recorded. Then, 
the value of the total losses (as a percentage of the number 
of photons in the input signal pulse) can be expressed as 
follows : 

C,{f^,L) = (1 - n^ff(f^, L)) ■ 100%. (48) 

Fig. [6] shows the losses associated with leakage C (blue 
curves, dotted lines) and the total losses of photons Cc, 
(red curves, full lines) as a function of the duration of 
the writing for a given medium length, and three different 
values of r. 

First, let us note that these curves are not monotonous 
and exhibit one or several minima. This means that for 
a given length L one can find the pulse duration which 
is recorded in the atomic medium with minimal losses. 
Significant difference between the curves corresponding to 
leakage and to total losses come from the role of the up- 
per level in the interaction of such pulses with the atomic 
medium. However, in the region of minimum losses, the 
distance between the two curves is small. The optimiza- 
tion of the memory based on leakage allows to define a 
range of values , for which an efficient writing is ex- 
pected. However, the curve giving the total losses is a more 
precise tool to determine the optimum ratio between 
and L. As can be seen from the plots, above some value of 
1^ the two curves coincide, i.e. all the system losses are 
associated with leakage only and level of |3) is not popu- 
lated at the end of the process. The larger detuning r the 
smaller the value for which this happens. 

We can also follow the dependence of the losses on the 
length L for a given value of as shown in Fig.]?] These 
curves have a monotonous variation, and show that the 
efficiency increases with increasing medium length. When 
the pulse duration increases, the curves of C and get 
closer to each other, in agreement with the analysis of Fig. 
[6) A specific feature of these curves is their saturation 
for large values of L. The presence of a plateau on the 
leakage plots is actually due to an approximation made in 
the model; we have neglected the time intervals associated 
with the propagation of the pulse wavefronts inside the 
medium, and we assumed i = is the time at which the 
wavefront reaches the the cell output, while t — is 
the time at which the end part of the pulse arrives at the 
entrance of the cell. This means that we always have some 
leakage in the initial time independently of the length of 
the medium. 

The difference between the levels of the plateau for £ 
and Cc characterizes the losses due to the non zero popu- 
lation of level |3). One can see that this value is constant 



for large enough length L. This value saturates because of 
the depletion of the signal field, so that further increase of 
the medium length cannot change the populations N2 and 
N3. The difference between C and Cc depends strongly on 
the pulse duration. In particular if the pulse is too short, 
many atoms are left in the upper state. It can also be 
seen in Figs. [7]: that when the detuning increases, the 
saturation occurs at larger values L, and that a high effi- 
ciency can be reached as well if the atomic medium is long 
enough. Various optimizations procedures have been pro- 
posed in the limit when the excited state can be adiabat- 
ically eliminated [21,22 . In particular, the optimization 
used in reference ^22j allows to get the maximum avail- 
able efficiency for long enough durations of the writing 
process, but breaks down when the duration of a writ- 
ing process gets smaller than the excited state decay 
time divided by the optical depth of the medium d. 
In Ref. [241 the numerical optimization procedure relying 
on the shaping of the driving pulse was extended to the 
non adiabatic case, which allowed to reach better storage 
efficiency for short pulses. The latter technique was de- 
veloped for the resonant case, yielding optimal memory 
efficiencies that are very close to the ones presented here. 
The applicability of our optimization method to various 
detunings shows that such a memory can be also very ef- 
ficient in the off-resonant regime, bringing more flexibility 
for experimental realizations. 

4.2 Validity limits of resonant and Raman 
approximations 

The solutions that we have presented are valid for detun- 
ings ranging from zero to large values that correspond to 
the case of Raman interaction, where the system is effec- 
tively reduced to a two-level system. General formulas cov- 
ering the full range of detunings allow a comparison with 
the limit cases of resonant and Raman interactions. We 
can identify the largest detuning for which the resonant 
approximation is still valid, yielding the same storage effi- 
ciency. On the other hand, we can estimate for which value 
of Zi a simplified Raman description can be used without 
yielding appreciable errors in the memory efficiency. 

In Fig. [2]d one can see a significant distortion of the 
temporal profile of the atomic ground state coherence at 
the input of the medium for r = 0.1 as compared to the 
resonant case. However this detuning does not affect the 
writing efficiency in a significant way. Figure shows the 
dependence of the total losses on time for a given length in 
two cases: for exact resonance (blue curve, dotted line) and 
with a detuning r = 0.1 (red curve, full line). The curves 
coincide to within 1.5% over the full range of . Thus, 
despite the local differences in the field-atom interaction in 
these two cases, the presence of a small detuning does not 
actually change the properties of memory cell as a whole. 
However when the detuning increases, the difference be- 
tween the curves increases (at r = 0.2 it reaches 4.5%, see 
Fig. m)), but in the range of interest for T^, that is the 
one that allows minimization of the losses, the curves are 
still close to each other (they agree within 1.5%). Further 
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increase of the detuning distorts the profile of the losses 
further, and the value of that provides minimum losses 
is shifted (see Fig. [8]:,d). 

Let us now turn to the case of large detuning. One can 
see in Fig. [9] that for r = 3 the profile of the total losses 
calculated with the general formulas (33)-(|34|) coincides 



well with the profile for the same quantity, calculated in 
the Raman approximation with ?' 3> 1 (the difference is 
about 2.5%, and it is less than 1% at the minimum). For 
a smaller detuning (for r = 2) the difference between the 
curves increases up to 7% (about 3% at the minimum): 
the calculation made in the Raman case underestimates 
the minimum losses and shifts toward lower values of . 
Thus, we can conclude that for ?■ = 3 and higher the Ra- 
man approach is applicable with good accuracy, but the 
general solutions should used for lower values of the de- 
tuning. 

For large enough detunings where the adiabatic limit 
is valid as well as our model, we can compare the results 
further. Our model predicts a storage efficiency below the 
maximal available efficiency reached by the method pro- 
posed in reference [22]. In the third plot of Fig. [6] the effi- 
ciency is 65% for Td-y = = 30 and d = 2400. For 
the same parameters the adiabatic optimization method 
converges to the maximal efficiency which is close to 100%. 
This is due to the fact that we do not elaborate shaping of 
a control pulse, used in reference 22^. Thus our method, 
even without control pulse shaping, is quite powerful in 
a non-adiabatic limit. Otherwise, the numerical adiabatic 
optimization of the control pulse profile should be used to 
reach the maximal storage efficiency for the long pulses. 




Fig. 9. Comparison of the general calculations (blue curves, 
dotted lines) and calculations in the Raman limit (red curves, 
full lines). Writing process : relative total losses (in percent 
of the input field intensity) at the output of the medium as a 
function of for L = 10 for (a)r = 2 and (b) r = 3. 



5 Discussion of the read-out process 

In Ref. [25j the optimization of the retrieval efficiency was 
based on the choice of the pulse duration providing the 
minimum leakage for a given medium length. Here we will 
look for a minimization of the total losses. Moreover, as 
was shown in Ref. [^[^ , backward retrieval provides sig- 
nificantly larger efficiency than forward retrieval, therefore 



we will study the cases of forward and backward retrieval 
and compare them. Using the result of Fig. [6^, we choose 
a value for input signal duration of =5.5 that provide 
minimum total losses for the writing process. 

For this value of we calculate the intensity of the 
retrieval field (normalized to the intensity of the input 
signal) as a function of the reading time (see Fig. 10 1). 
The plot obtained from optimization based on leakage |25j 
is shown in Fig. [T0|3, in order to compare the results. The 
retrieval efficiency is defined by: 



, , f' |a^(i,z,q)|2di 
£(q) = ^^^^J ^ ^' X 100% 

Jo \ain{t,q)\'^dt 



(49) 



In the case of optimization based on total losses, we 
find that the retrieval efficiency is equal to 88% at = 
2T^ . In view of the writing efficiency {Cc — 11.2%), we 
see that using a reading time = 2f^ we can restore 
almost all the information written in the medium. When 
the leakage-based optimization is used, we find £ — 84% 
at — 37'^. This comparison clearly demonstrates the 
benefits of the optimization based on the total losses. 

As for the temporal profiles of retrieval field, it is ob- 
vious that in both cases they are very different from the 
input signal profile, which is a usual result in most mem- 
ory processes. 

Let us now consider the result of the forward and back- 
ward retrieval for the non- resonant case (see Fig. 11 for 
r = 0.5). 

The efficiency is much lower for forward retrieval than 
for backward retrieval one: at = 10 (i.e. « 3.37'^) 
the efficiency of forward process is £ = 58%, while the 
efficiency of backward process is £ = 85.6%. Moreover, 
in the latter case less than 2% of the available photons 
remain in atomic ensemble (since the writing losses were 
12.6%). Comparison with the resonant case shows that the 
presence of detuning slows down the read-out. 

Finally, we can show that even for a large detuning, 
the total efficiency can be large, at the condition that the 



medium length L is large enough. In Fig. 12 the red (oscil- 
lating) curve corresponds to the calculated read-out signal 
intensity for backward retrieval for L = 100, r = 2, and 
for a pulse duration of the input field = Att, which 
ensures a minimum in the total writing losses (writing ef- 
ficiency is n'^-f^ — 95.6%). The oscillation period is equal 
to 47rr, similar to the modulation observed in Fig. [2] If 
we calculate the intensity of the retrieved signal in the 
Raman approximation, the overall shape of the curve re- 
mains the same, but the oscillations disappear (blue curve. 



dotted line in Fig. 12 1. The retrieval efficiency calculated 



with these curves differ by 1.4% (78.6% for the exact cal- 
culation and 77.2% for the calculation in Raman approx- 
imation), so that the calculation in Raman limit can be 
a quite good estimation for memory efficiency at r = 2. 
Note, that here like for the curves in Fig. [2] the magnitude 
of the oscillations will decrease with increasing detuning. 

In addition to high speed operation, our model includes 
transverse coordinates, as can be seen from Eq. ([2]), and 
thus allows the treatment of a variety of multimode fields. 
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Fig. 12. Reading process : field intensity \a^{i,L)\^ at the 
output of the medium for L = 100, TV ~ 47r and r — 2 for 
backward retrieval calculated with exact formulas (red curve, 
full line) and within the Raman approximation (blue curve, 
dotted line). 

in the same way as in Ref. [55] . The storage of a multimode 
signal field is performed using a single mode signal field, 
and the information is written in the atomic ensemble 
as a hologram. The quantum hologram process together 
with a Raman memory was proposed for the first time by 
Sokolov et. al. in Ref. [2D] and it was shown that the mem- 
ory capacity is limited by diffraction, but with the limita- 
tion depending on the direction of readout. Under forward 
readout the maximum number of the modes or picksels N 
which can be stored in principle is given by the square 
of the Fresnel number N ~ F^, where Fjy = S/{XL). 
This expresses the condition that the output pixel size 
should not exceed the transverse size of the memory cell. 
This rather loose limitation comes from the fact that most 
diffraction effects are compensated between writing and 
readout with similar geometries. A similar compensation 
does not take place for the backward readout and as a 
result the number of stored modes TV in this case can not 
exceed the Fresnel number. 



cally not involved in the process, all three levels are popu- 
lated, and this leads ultimately to additional undesirable 
losses. We have shown that it is possible to choose the 
parameters to make these losses very small. By control- 
ling the Rabi oscillation one can preferentially populate 
the lower state |2) and depopulate state |3). When the 
frequency detuning increases, the upper state population 
naturally decreases, which contributes to the reduction of 
losses. Moreover, for large frequency detunings, the in- 
teraction between the atomic medium and the fields is 
weaker, and a larger medium depth is necessary to reduce 
the leakage. 

An important question is how to optimize the memory 
process to obtain the highest possible quantum efficiency. 
In the literature two approaches are mainly discussed. The 
first one is based on the choice of the optimal time-shape 
for the signal pulse, using for this the eigenf unctions of the 
integral operators specific of the considered memory [21] . 
It was demonstrated that if the functions possess some 
symmetry, the efficiency can be improved. In our case with 
very short pulses, there is no such symmetry and this op- 
timization turns out to be impossible. Another optimiza- 
tion method was proposed by Gorshkov et. al. in Ref. [22], 
based on a study of the driving pulse shape. However, 
the search procedure for the optimal shape in this arti- 
cle is based on the adiabatic or Raman approximations, 
which are not generally applicable in the case of very short 
pulses. Non-adiabatic pulse shape optimization was also 
developed for the resonant case [24] extending the validity 
of the method. In our model, the optimization is based on 
the minimization of the losses, which come from the leak- 
age and from the population of the upper atomic level, 
and allows identifying the optimal combination of signal 
pulse duration and optical depth. We demonstrate that 
high memory efficiencies can be achieved by this method 
for very short pulses whatever the value of the detuning. 



6 Conclusion 

In this article, we have presented a model for a high speed 
quantum memory based on a three-level medium in the 
yl-configuration. This model relies on a full calculation of 
the fields and of the atomic coherences as a function of 
time and space for arbitrary frequency one-photon detun- 
ing, while keeping the two-photon resonance. It allows to 
examine all the situations between interactions resonant 
with the excited state and non resonant interaction, which 
corresponds to a Raman transition. Our model allows to 
identify the conditions in which the interaction can be 
treated as a Raman transition. This corresponds to a nor- 
malized detuning r — zi/(2i?) > 2. For larger detunings, 
we have shown that the Raman model gives accurate re- 
sults. 

In the near-resonant case, the interaction of the three- 
level atomic system with the weak signal field and the 
strong driving field turns out to be more complicated than 
in the Raman case. As a matter of fact in contrast to the 
Raman process, where the upper atomic state is practi- 
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A General solutions for the main equations 



Eqs. (|26|-(28) can be solved in the general form by using 



the Laplace domain. Some details of the formal proce- 
dures are discussed in App. [B] and now we start with the 
solutions themselves in the explicit form. For our analysis 
in this article we do not need full information about the 
solutions, nevertheless one can find it below. Here and ev- 
erywhere in article under consideration of the solution we 
use dimensionless co-ordinates t and z given by 



t = m, 



2g^N 



(50) 
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Let us write the general solutions in the form under the 
arbitrary initial and boundary conditions 

a(t,z;q)= / di' ain{i - i' ; q)Gaaii' , z) 



1 

'2pi 



di'b{0,z- z';q)Gbaii,z') 
dz'c(0,z-z';q)Gca(t,2'), (51) 



b{t,z]q)= -p / dt'ain{t-t';q)Gabit',z) 
Jo 

+ 1 f d5'b{0,S~S';q)GbbilS') 



+ i / dS'c{0,Z-S';q)Gcb{i,S'), (52) 



c(t, i;q)= p / dt'din{t-t':q)Gac{t',z) 
Jo 

1 



+ - / di'b{0,2~S';q)Gbcit,i') 
^ Jo 

+ i / dz'c(0,z-z';q-)Gcc(t,5'), (53) 
^ Jo 

where kernels Gik{t,z) are bilinear combinations of the 
expressions depending on the n-th Bessel functions of the 
first kind denoted by J„ 



fit,i;r)=6{t)-e 

X Ji (^^(l + r)^^ 0{i) 



+ /(I + r)5 

V 4! 



(54) 



fi{i,z; 



foii,z;r\ 



(VT+72 + r)i U{1 + r)i 



xJi (^^(l + r)tzj e(i), 



(55) 



xJo V(l + ^)^^^(^) 



The kernels in Eq. (51 1 read 



Gaa{i,S)^[.f{r)^ri-r)] (t, 5), 
Gb,(t, f) = [/o(r)*/*(-r)](i,S), 

G,a(i,^) = ^[/o(0 * r{-r)]{i,~z) 



(59) 



We denote the t convolution of two arbitrary functions 
X{i,z;r) and Y{i,z;r) as 

[X{r)*Y*{-r)]{irz) = 

rt 

di'X{i- i', z; r)Y*{i', 5; -r). (60) 



The parameter r is the dimensionless frequency detuning 

r = AI{2n). 

The other kernels read 

Gab{t,i) = [h{r) * f*{-r)]{i,i), (61) 
Gbb{l z) - 2 5{z)F,{t) + [/i(r) * /r(-r)](t, z), (62) 

G,6(f, 5) = 2 <5(5)F2(i) + ^[/i(r) * /o*(-r)](i, z) 



(63) 



and 

Gac{t, z) 



1 + r 
2 

1 - r 



[/o(r)*r(-r)](f,z) 



[/('^)*/o (-01(^,5) 



(64) 



Gbc(i,^) = 2(5(5)F2(t) 
1 - r 



1 + r 



[hir)*f^{-r)]{irz) 



[/o(r)*A*(-r)](t,5) 



(65) 



GUi, S) = 2 S{z)Fs{i) - ^^^[/i(r) * /*(-r)](f, S) 



(1- 



-[/(r)*/r(-r)](t,z) 
[/o('-)*/o*(-r)](?,5)- 



(66) 



In these formulas the time dependent factors F{i) read 



Flit) 
F2{i) 



cos I V -L + 7" t 
ir 

: sm 



1 



: sm 



v/l + r2 i 



cos (^VTTr^ 



sin ^VTTr2 tj 



-iri 



(67) 
(68) 

(69) 



(56) B Solutions using the Laplace transform 



(57) 
(58) 



In order to solve Eqs. (|26[)-(|28[), we rewrite them in the 
Laplace domain 



dz 



ds{z; q) = -gVNcs{z; q) 



(70) 



-c(z,0; q) + (s + iA)Cs{z]q) = gVNas{z:q) + nbs{z;q), 

(71) 

~b{z, 0; q) + sbsiz; q) ^ -/2c,(z; q). (72) 

From this we can write the closed differential equation for 
the field amplitude as{z] q) in the form 

dds{z\ q) 
dz 



rsas{z;q) ~ gVNas{z;q). (73) 
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Here the coefficient Fg determines a rate of escape of the 
amplitude along the z-axis and reads 



(74) 



2 \s + ifiQ s~ivf}J^ 
where the following notations are introduced 

f} = n^/l + r\ = 1 + r, 1/ = 1 - r. (75) 



The inhomogeneous term on the right in Eq. ( 73 1 is deter- 



mined by the initial conditions for the medium state and 
given by 



as(z;q) 



1 



s{s + iA) + Q^ 



i76(0,z;q) + sc(0,z;q) .(76) 



A solution of Eq. ( 76 1 reads 



a,{z;q) = a,iO;q)e-'^^^'- -gVN J dz'a,(z'; q)e-^=(^-^'). 



(77) 



From Eqs. (71)-(72| one can obtain 
s 



Cs{z;q) gVN 



1 



s{s + iA) + 



as{z;q) +as{z;q), (78) 



bs{z;q)^- h(Q,z\q) - Qcs{z;q) 



(79) 



After the inverse Laplace transformation in Eqs. (77)- ([79 
one can obtain all the solutions in the form 



a{t,z;q)^ / dt' ain{t - t';q) D{t' ,z) 



(80) 



-gVN dt' I dz'a{t-t',z- z';q) D{t',z'), 
Jo Jo 

c{t,z;q) = 

gVN dt'F3{t-t')a{t',z;q) + ait,z;q), (81) 



b{t,z:q) = 



gVN / dt'F2{t ~ t') d{t' ,z;q) + f3{z,t;q), (82) 



where the terms a{t, z; q) and j3{z, t; q) are expressed via 
the initial conditions for the medium coherences in the 
form 

ait, z; q) = F^it) c(0, z; q) + F2{t) 6(0, z; q), (83) 
/3(z, t; q) = -F^it) c(0, z; q) + F,{t) S(0, z; q), (84) 

and the kernel D{t, z) is simply proportional to Gaa (see 
(56)) and is given by 



D{t,z) = nGaa{i,i). (85) 
After some simple transformations one can obtain the 



solution in the form (51)-(53). At the same time for the 



References 

1. A. I. Lvovsky, B. C. Sanders, and W. Tittel, Nature Pho- 
tonics 3 (2009) 706 

2. C. Simon, M. Afzelius, J. Appel, A. Boyer de la Giroday, 
S. J. Dewhurst, N. Gisin, C. Y. Hu, C. Y. Jelezko, S. KroU, 
J. H. MiiUer, Eur. Phys. J. D 58 (2010) 1 

3. K. S. Choi, H. Deng, J. Laurat, and H. J. Kimble, Nature 
452 (2008) 67 

4. K. Honda, D. Akamatsu, M. Arikawa, Y. Yokoi, K. Akiba, 
S. Nagatsuka, T. Tanimura, A. Furusawa, and M. Kozuma, 
Phys. Rev. Lett. 100 (2008) 093601 

5. J. Appel, E. Figueroa, D. Korystov, M. Lobino, and A. I. 
Lvovsky, Phys. Rev. Lett. 100 (2008) 093602 

6. C. Clausen, I. Usmani, F. Bussieres, N. Sangouard, M. 
Afzelius, H. de Riedmatten, and N. Gisin, Nature 469 
(2011) 508 

7. E. Saglamyurek, N. Sinclair, J. Jin, J. A. Slater, D. Oblak, 
F. Bussieres, M. George, R. Ricken, W. Sohler, and W. 
Tittel, Nature 469 (2011) 512 

8. J. Cviklinski, J. Ortalo, J. Laurat, A. Bramati, M. Pinard, 

E. Giacobino, Phys. Rev. Lett. 101, (2008) 133601 

9. Hugues de Riedmatten, Mikael Afzelius, Matthias Staudt, 
Christoph Simon, Nicolas Gisin, Nature 456, (2008) 773 

10. Hedges, M., Longdell, J., Li, Y. & Sellars, M., Nature 465, 
(2010) 1052 

11. M. Hosseini, G. Campbell, B. M. Sparkes, P. K. Lam and 
B. C. Buchler, Nature Physics 7, (2011) 794-798 

12. K. F. Reim, P. Michelberger, K. C. Lee, J. Nunn, N. K. 
Langford, and I. A. Walmsley, Phys. Rev. Lett. 107, (2011) 
053603 

13. M. Fleischhauer and M.D. Lukin Phys. Rev. Lett. 84, 
(2000) 5094 

14. G. Hetet, Longdell, J. J. Alexander, P. K. Lam, and M. J. 
Sellars, Phys. Rev. Lett. 100, (2008) 023601 

15. S. A. Moiseev and S. KroU, Phys. Rev. Lett. 87, (2001) 
173601 

16. M. Afzelius, C. Simon, H. de Riedmatten, and N. Gisin, 
Phys. Rev. A 79, (2009) 052329 

17. W. Tittel, M. Afzelius, T. Chaneliere, R. L. Cone, S. KroU, 
S. A. Moiseev, and M. Sellars, Laser & Photonics Reviews 
4, (2010) 244 

18. M. Hosseini, B. M. Sparkes, G. Campbell, P. K. Lam, and 
B. C. Buchler, Nature Communications 2 174 

19. D. V. Vasilyev, L V. Sokolov, and E. S. Polzik. Quantum 
memory for images: a quantum hologram. Phys. Rev. A 
77, (2008) 020302(R) 

20. D. V. Vasilyev, L V. Sokolov, and E. S. Polzik. Quantum 
volume hologram. Phys. Rev. A 81, (2010) 020302(R) 

21. J. Nunn, L A. Walmsley, M. G. Raymer, K. Surmacz, 

F. C. Waldermann, Z. Wang, and D. Jaksch, Phys. Rev. 
A 75,(2007) 011401 

22. A. Gorshkov, A. Andre, M. D. Lukin, and A. S. S0rensen, 
Phys. Rev. A 76 (2007) 033805 

23. K. F. Reim, J. Nunn, V. O. Lorenz, B. J. Sussman, K. C. 
Lee, N. K. Langford, D. Jaksch and L A. Walmsley Nature 
Photonics 4, (2010) 219 

24. Alexey V. Gorshkov, Tommaso Calarco, Mikhail D. Lukin, 
and Anders S. S0rensen, Phys. Rev. A 77, (2008) 043806 

25. T. Golubeva, Yu. Golubev, O. Mishina, A. Bramati, J. 
Laurat, and E. Giacobino, Phys. Rev. A 83, (2011) 053810 

26. Mikhail L Kolobov, Rev. Mod. Phys. 71, (1999) 1539 



numerical computation it is possible to use Eqs. (80l-(82) 



T. Golubeva et.al.: High speed yl-type memory with arbitrary frequency detuning 



11 



\b>'m\' 




\h^(I.O)\' 0.6 



Fig. 2. Normahzed coherence at z = as a function of time for (a) r = 0,(b) r = 0.1, (c) r — 0.2, (d) r = 1, (e) r = 2, (f) 
r = 10. 






Fig. 3. Distribution of the coherence \b^{t,z)\^ in time and space for (a)r- — 0, (b) r — 1 and (c) r — 2. 
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Fig. 4. Normalized distributions of the atomic coherence inside the medium at time 
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t = TT for (a)r = 0, (b) r = 0.5, (c)r = 2. 
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Fig. 5. Fraction of signal photons (normahzed to the energy of input signal pulse) that have been converted to atomic coherence 
during writing as a function of the detuning parameter r for L = 10 and (a) — 7r/4, (b) — 7r/2, (c) = tt, (d) 
f ^ = 2-K. 
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Fig. 6. Writing process : relative losses associated with leakage (blue curves, dotted lines) and relative total losses (red curves, 
full lines) (in percent of the input field intensity) at the output of the medium as a function of Tw for L — 10 for (a)r — 0, (b) 
r = 0.5, (c)r = 2. 
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Fig. 7. Writing process : relative losses due to leakage (blue curves, dotted lines) and relative total losses (red curves, full lines) 
(in percent of the input field intensity) at the output of the medium as a function of L and (a) r — 0, (b) r = 0.5 and (c) r = 2 
for Tw = 7r/2 (first row) and Tw ~ tt (second row). 
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Fig. 8. Comparison of relative total losses (in percent of the input field intensity) at the output of the medium as a function of 
Tw for resonant (blue curves, dotted lines) and detuned (red curves, full lines) cases; L = 10, (a)r = 0.1, (b) r = 0.2, (c)r — 0.3, 
(d)r = 0.4. 
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Fig. 11. Reading process : field intensity |a^(f, L)p at the output of the medium for L — 10, Tw ~ 3 and r = 0.5 for (a) 
forward and (b) backward propagating retrieval. 



